Molecular weight effects on chain pull-out fracture of reinforced polymeric interfaces 
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Using Brownian dynamics, we simulate the fracture of polymer interfaces reinforced by diblock 
connector chains. We find that for short chains the interface fracture toughness depends linearly 
on the degree of polymerization A'^ of the connector chains, while for longer chains the dependence 
becomes N^^^ . Based on the geometry of initial chain configuration, we propose a scaling argument 
that accounts for both short and long chain limits and crossover between them. 
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Creating advanced materials often means mixing dif- 
ferent homo-polymers to produce systems with desired 
combined properties. However, most polymer blends are 
immiscible: They form macroscopically phase-separated 
mixtures with only interfacial van der Waals forces keep- 
ing domains of different phases together. The fracture 
toughness of such blends is limited by one of the inter- 
faces: In ideal defectless conditions, it is equal to the 
work of adhesion W — + — Jab between two homo- 
polymer phases. Here, ja and are respectively the 
surface energy of homo-polymer A and homo-polymer 
B and jab is the interfacial free energy. Reinforcement 
of these weak polymeric interfaces is often achieved by 
the addition of A-B diblock copolymers which compati- 
bilize the blend and strengthen the interface The 
strengthening can be attributed to the miscibility of each 
block with one of the homopolymers. This causes the 
block copolymer to expand and entangle with homo- 
polymer phases on either side of the interface. The inter- 
facial tension (energy) is reduced, the interfacial width 
is increased, and the adhesion thereby improved. Frac- 
ture toughness and failure mechanisms of such reinforced 
polymeric interfaces have been investigated extensively 
by, for example, experiments on different incompatible 
systems of polymer glasses [^-p^ and cross-linked net- 
works (elastomers) [ pl] , p^ . 

Several theoretical models have been proposed 13-1^] 
to explain the reinforcing effect of connector chains in 
both elastomers and glassy polymers. A "failure mecha- 
nism map" has been developed 15,1^,|I^ which relates 
the mechanism of interface failure to the polymerization 
index N , surface density a of connector chains, and the 
time scale on which the deformation occurs. According 
to the failure map there are three major mechanisms: (i) 
chain scission, which happens whenever the stress along 
the connector chain becomes larger than the strength of 
the covalent bond between segments of the chain, (ii) 
pull-out of the connector chain as a result of disentan- 
glement from homo-polymer phase, and (iii) failure by 
craze formation, followed by chain scission or chain pull- 
out, which take place when a large stress is transfered to 
the bulk of the homo-polymer phases. 



In this paper we direct our attention to the case where 
interface failure is due to chain pull-out, and specifi- 
cally focus on the effect of the polymerization index N of 
the connector chain. We consider a "mushroom regime" 
p"8[-p3[ where connector chains are grafted with low sur- 
face density, cr <C 1/NP {N is the number of monomers 
per grafted chain and I is the monomer size). In this 
low-density regime, the equilibrium shape of the connec- 
tor chain in the homo-polymer phase is mushroom or 
plume-like. Furthermore, mutual entanglement between 
different connectors is negligible. Hence we need con- 
sider only the behavior of a single chain. Experimental 
data on the dependence of G on iV is scattered, typically 
assuming a form 



G- A", 



(1) 



where estimates give variously 1 < a < 2 pUlCj]. Both 
linear and quadratic dependences of G on A have been 
predicted, utilizing a tube picture |^^. Different con- 
stitutive equations have been proposed which relate the 
local stress in planar cohesive zone near the crack tip to 
various phenomena, such as, the penetration depth of the 
chains, their surface density, and the pulling rate [^-^. 
These studies have been focused mainly on the pull-out 
fracture in tensile mode, when the applied force is nor- 
mal to the interface. An alternative mode of interface 
failure is shear fracture, or the resistance of the inter- 
face against slip. In fracture mechanics of bulk materials 
generally only the tensile or opening mode is important, 
since cracks normally travel in a direction that maxi- 
mizes the opening mode. However, for the interface be- 
tween two different materials, the situation can be more 
complex: The crack is often constrained to follow the in- 
terface, giving rise to the possibility of crack propagation 
involving a combination of tensile and shear modes 
Herein we study the dynamics of chain pull-out fracture 
in both tensile and shear mode separately and examine 
the dependence of G on A. 

We use a Brownian dynamics method, due to Picket, 
Jasnow, and Balazs to simulate the pull-out of a sin- 
gle connector chain of length A from a two-dimensional 
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homo-polymer phase. To quantify the interface tough- 
ness for different chain length, we calculate the work 
that is required to pull out the chain with constant ve- 
locity vq. The homo-polymer phase is modeled by a 
two-dimensional semi-infinite square lattice of obstacles. 
Each obstacle represents an entanglement or cross linked 
point, depending on the glassy or elastomeric structure 
of the homo-polymer phase, and provides lateral con- 
straint on the movement of the connector chain. The 
connector is represented by a freely jointed chain with N 
links wherein there is no self-interaction. Hence individ- 
ual monomers can freely pass over each other, provided 
links have a constant length. The initial configuration 
is created by putting the first monomer at the bound- 
ary and then continuing the chain as a random walk in 
dimension c? = 2, with a reflecting boundary at the in- 
terface. The random walk is restricted by requiring, as 
we shall see below, that monomers are repelled by the 
obstacles. The pull-out dynamics is simulated by pulling 
the chain by the first monomer at a constant velocity in 
either tensile mode (perpendicular to the interface) or 
shear mode (parallel to the interface). This is done con- 
veniently by moving the obstacle matrix at velocity — vq 
while keeping the first monomer fixed. That is, if ri is 
the position of the ith monomer, dri/dt = 0. 

The movement of the other i = 2, . . . , iV monomers is 
governed by the following over-damped Langevin equa- 
tion. 

1^ - vo^ = Tj(ri+i - ri) - Tj_i(ri - ri_i) + Fi + rji, 

(2) 

where fn+i = tn- Here is the amplitude of ten- 
sion in the segment connecting monomers i and i -I- 1, 
r]i is the random Gaussian noise representing the effect 
of thermal fluctuations on each monomer, and is the 
viscous friction coefficient per monomer. The strength of 
the random noise is related to the monomeric friction 
by the fluctuation-dissipation theorem {rii{t)rij{t )) — 
2vkBT5ij5{t — t ), where ks is Boltzmann's constant and 
T is temperature. Fi(Ari) is a short range monomer- 
obstacle repulsive force, where Arj is the distance be- 
tween the ith monomer and the closest obstacle. For 
I Aril < J'c, there is hard-core repulsion. For Tc < |Ari| < 
2rc, there is soft-core repulsion, where the force obeys 



(4) 



Fi = o-Ari 
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|Ari|2 
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(3) 



The parameter a controls the strength of the force [gSj . 
For larger separations, |Ari| > 2rc, there is no force act- 
ing on the monomer. In our simulations we have used 
Tc = 1/2, where I is the link length of the connector 
chain. 

To solve Eq. (||) , the tensions in each segment of the 
connector chain have to be determined. This is done by 
enforcing the constant segmental length constraint: 



for 1 < i < iV — 1. This results in a tridiagonal matrix 
equation for the r^'s [E5|: 



^ ' = n+i(ri+2 - ri+i)(ri+i - n) - 2nf + 

+ Ti-i(ri - ri_i)(ri+i - ri) + 
+ (Fi+i + r/i+i - Fi - 77i)(ri+i - n). (5) 



Ideally, by applying the constant segmental length con- 
straint, the left-hand side of Eq. (4) should be zero. How- 
ever, after several updates of the simulation, as a result 
of accumulation of roundoff errors, the distances between 
some of the adjacent monomers differ slightly from I. 
To improve the stability of the algorithm, a correction, 
restoring the original segment length, is introduced: 



d(ri+i - ri)^ 
2dt 



(|ri+i-ri|^0 
dt 



(6) 



Eq. (g) is solved using a fourth-order Runga-Kutta algo- 
rithm. 

We consider reduced units in which i/, wq, cr, /csT, and 
D, the distance between obstacles, are set equal to unity. 
The distance between two adjacent monomer in connec- 
tor chain I is chosen to be 0.4. The time step of the 
simulation, typically lO"** to 10~^, is adjusted so that 
the average difference in segment length from I in each 
run is less than 0.1 percent. The results were averaged 
over 20 — 100 independent realizations of the initial con- 
ditions. 



Figure 1 (a) 
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Figure 1{b) 
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FIG. 1. Snapshot of a simulation of fracture for a chain of 
A?" = 700 monomers: (a) tensile and (b) shear mode. 



tensile and shear modes. Naturally, the fracture tough- 
ness grows as the degree of polymerization of the con- 
nector chain increases. This growth is more significant 
for larger N due to the fact that longer connectors can 
penetrate well beyond the neighborhood of the interface. 
They entangle efficiently with the bulk polymers of the 
compatible phase. It also shows that the effect of rein- 
forcement on the interface is higher in tensile mode, by 
a factor of about five for the present model ||2^. One 
can notice the existence of two scaling regimes for G vs. 
N, corresponding to short connectors, N < N^, and long 
connectors, N > Nc where Nc is the crossover length. 

Figure 3 shows the power-law dependence of G on N. 
For short connectors G scales linearly with N, G ^ N, 
while for long connectors it scales as G ~ ]\[^/^ . 



Figure 3 



G~N" 



Fig. 1 shows a snapshot of the simulation for both ten- 
sile and shear pull out. We quantify the fracture tough- 
ness of the interface by determining the work G required 
to remove the connector chain completely from the ob- 
stacle lattice. The work is G = / Pdt, where the power 
P is obtained via P — vqt, vq is the constant pulling 
velocity, and r is the instantaneous tension in the seg- 
ment crossing the boundary of the obstacle lattice in the 
pulling direction at each time step. 



Figure 2 
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FIG. 3. Log-log plot of G vs. A'^ for long connectors 
N = 150 to 900. (•) and (o) represent tensile and shear 
modes respectively. Data is consistent with G ~ N^^^ (solid 
line) . 
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FIG. 2. Ln-Ln plot of G vs. A''. (•) represents tensile and 
(o) indicates shear mode. Solid lines show that short-chain 
data (TV < 100) are consistent with a = 1, while long-chain 
data (TV > 100) are consistent with a = 3/2. 

Fig. 2 shows the fracture toughness G as a function of 
the polymerization index N of the connector chain for 



These regimes can be understood as follows. For small 
TV, the chain is entangled with one or two layers of ob- 



stacles, even if the free radius of gyration Rg 
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is of order or smaller than the obstacle lattice spacing 
D. Since the the initial configuration is repelled by the 
obstacles, the available phase space for a chain is quite 
restricted. For the first few chain segments, the reflect- 
ing boundary condition and the the area taken up by the 
obstacles, results in the chains going straight into the lat- 
tice. After passing the first row or two of obstacles, the 
chain, if it is long enough, usually bends, with tight loops 
around obstacle cores being highly unprobable. When 
the chain is pulled out, the dominant deterministic part 
of the tension equals the total viscous drag force. This 
force is roughly proportional to the number of monomers 
Nmove simultaneously in motion. Since pulling tension 
cannot propagate through a loop unless it is tightened 
around an obstacle core, for short chains Nmove is of or- 
der of the number of monomers stretched in one obstacle 
lattice spacing D, that is D/l. Hence for small TV, the 
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total work to pull the chain out can be estimated as 

Tj rNl/vo 

G « vo"^!^— / dt^ VQi^ND, (7) 
' Jo 

which is consistent with the linear scaling regime ob- 
served in simulations for short connectors. 

For large N, the penetration depth increases. Indeed, 
for sufficiently large N, the chain center of mass is lo- 
cated at a distance of order of the radius of gyration Rg 
from interface. Then, the number of monomers that are 
simultaneously in motion becomes proportional to this 
distance, Nmove ^ Rg ^ N^/"^, which we have observed 
directly. Consequently, the pull-out fracture energy be- 
comes 

fNl/vo 

G^VO Fdragdt « VoI^N'^^H. (8) 

Jo 

This is roughly analogous to what happens when one 
pulls on a wet garden hose of length N, left on the ground 
by a gardener who has performed a random walk among a 
grove of trees. Upon pulling, only the currently stretched 
segment starts to move, if friction in the self-intersections 
of the hose is negligible. If the hose configuration was 
created as a result of random walk, the length of such a 
segment scales as a typical size of such a walk, %/7V, where 
N is the total length of the hose. The hose configuration 
may include a few loops around the tree trunks, which, 
upon tightening, also start to move. The total length 
of these tightened loops has no important contribution: 
That length scales as a winding angle of a random walk 
of the length N, i.e., as IniV (see, for example, Ref. [pTf). 
As a result, the drag force is proportional to the length 
of the hose that is constantly in motion, ^/N. 

Although our numerical work is in two dimensions, 
these physical arguments also follow in three dimensions, 
and are applicable to experimental systems The 
small N regime crosses over to a large N regime when 
the radius of gyration Rg becomes of order of the obstacle 
spacing D. That gives Ncross [D/l]'^, and we expect 

G = Nf{N/N,ross). (9) 

where the crossover scaling function obeys f{x 0) = 
const, and f{x oo) = x^/^. 

As mentioned earlier, different forms of the power-low 
dependence of G on have been proposed in previous 
theoretical studies: Both G ^ N and G ~ A^^ have been 
predicted. The results from experiments are somewhat 
ambiguous. The linear dependence of G on A^ is 
predicted for very slow crack propagation velocities, and 
therefore corresponds to G — > Gq, where Go is the frac- 
ture toughness threshold or minimal energy required to 
break the interface when vq 0. On the other hand, 
G ~ A^^ is predicted for pull-out fractures where the 
crack propagation speed is high or at least larger than 
a critical velocity. In other words, these two predicted 



regimes differ by the pull-out velocity. The crossover be- 
tween the two scaling regimes G ~ A" and G ~ N^f^ 
obtained in our simulations takes place in the latter high 
velocity regime: The pulling velocity used in our sim- 
ulations Vq = ksT / {vD) corresponds to relatively fast 
pulling rate |^9|] , which in turn results in high crack prop- 
agation speeds. Our work is in the same regime, and 
indeed extends the numerical work and theoretical argu- 
ments of Picket, Jasnow, and Balazs |jl^, wherein this 
model was introduced. For large A, Picket et al. argued 
G ~ A^^ in contrast to our result of G N^^^. This 
was based on an analogy of the polymer to the motion 
of a rope in "block and tackle" pulleys. In that case, 
the average drag force to pull the connector chain out 
of the matrix is fdrag ~ vvqN. Hence, it follows that 
the number of monomers instantaneously in motion are 
Nmove ^ A. In fact, the correct analogy is not to pulleys, 
but to the motion of a garden hose lying on grass, as de- 
scribed above. Then fdrag ~ N^^^, giving rise to the 3/2 
exponent. Our numerical results, which are more exten- 
sive than those of the earlier work, support this picture. 

To conclude, we studied the chain pull-out fracture of a 
reinforced polymeric interface in tensile and shear modes 
using molecular dynamics algorithm. Our results confirm 
the nonlinear dependence of the fracture toughness of the 
interface on the length of the connector chain observed in 
experiments. We found that, depending on the length of 
connector chain, the fracture toughness of the interface 
shows different scaling dependences: For short chains, G 
scales linearly with A^, while for long connectors, we ob- 
served a crossover to a new scaling regime G ~ N^^^. 
Our results can be tested experimentally on elastomeric 
networks for long connector chains in the low coverage 
mushroom regime 
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